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Abstract. The spin 1 bilinear-biquadratic model H — Yluj) [ cos ' Sj + sin <p(Si ■ S 
i-^j ■ on square lattice in the region < <fi < ir/A is studied in a fermion representation with 

G \ a p-wave pairing BCS type mean-field theory. Our results show there may exist a 

non-trivial gapped spin liquid with time-reversal symmetry spontaneously breaking. 
This exotic state manifests its topological nature by forming chiral states at the edges. 
To show it more clear, we set up and solved a ribbon system. We got a gapless dis- 

T 1 

■ persion representing the edge modes beneath the bulk modes. The edge modes with 
£N1 . nonzero longitudinal momentum (k x =/= 0) convect in opposite directions at the two 

edges, which leads to a two-fold degeneracy. While the zero longitudinal momentum 

■ (k x = 0) modes turn out to be Majorana fermion states. The edge spin correlation 

functions are found to decay in a power law with the distance increasing. We also cal- 
culated the contribution of the edge modes to the specific heat and obtained a linear 

CnJ ■ law at low temperatures. 



PACS numbers: 75.10.Kt, 05.30.Fk, 74.20.Fg, 73.43.-f 
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1. Introduction 

In condensed matter physics, the Landau's theories of Fermi liquid and spontaneous 
symmetry breaking have been the basic principles that account for vast phenomena. For 
instance, the ground state of the two-dimensional spin S Heisenberg model on a square 
lattice possesses long range Neel order and the Goldstone modes due to spontaneous 
spin-rotation symmetry breaking. People have been trying to look for more exotic states 
in spin systems for a long time. The seminal concept of "resonating valence bond" 
(RVB) spin liquid state was first proposed by P. W. Anderson [I]. Then, mostly due to 
its implication of the mechanism of the high temperature superconductivity [2], this field 
has been flourishing for more than two decades. As a new type of quantum matter, the 
spin liquid state itself is intriguing since its properties have never been clarified before [3] . 
Various approaches have shown that quantum spin liquids may exist in two-dimensional 
(2D) S — 1/2 J1-J2 model and the Heisenberg model on the Kagome lattice. In these 
models, the quantum spin liquids are accessed (in principle) by appropriate frustrating 
interactions [4]. However, the nature of the quantum disordered ground state is still 
under debate. RVB spin liquid state obviously goes beyond the Landau's theories, in 
which the quasiparticles of Fermi liquid carry both spin and charge quantum numbers. 
Another exciting field in searching for quantum exotic states beyond the Landau's theory 
of spontaneous symmetry breaking is the quantum Hall (QH) and fractional quantum 
Hall (FQH) states. In these states, topological order plays an essential role |5J In 
the QH state, a quantized Hall conductance was measured due to the formation of the 
Landau levels of a 2D electron gas at low temperatures and in strong magnetic field. One 
of the key features of QH effect is the existence of chiral edge states around the system 
boundaries. Recently, involving both quantized spin Hall effect (TJ [8] and quantized 
anomalous Hall effect [9], the topological insulator with gapless edge states has been 
one of the hot issues. The fundamental links between the above two fields have attracted 
much attention [10]. In this work, we aim to contribute to these interesting topics. 

We study the spin 1 bilinear-biquadratic model on square lattice, of which the 
Hamiltonian is written as 



where Sj is a spin 1 operator. Its semiclassical version (S — > 00) on a bipartitie 
lattice exhibits four ordered phases that are exactly divided by four SU(3) symmetric 
points at the model parameters, = 0o — vr/4, ±7r/2, — 37r/4 [TTJ [12]. They are 
phases with ferromagnetic (FM), antiferromagnetic (AFM), ferroquadrupolar (FQ), 
and antiferroquadrupolar (AFQ) orders respectively. Whether these classically well- 
understood phases are stable in the quantum case or how the quantum model behaves 
constitues the interesting topics of current researches. In one dimension (ID), many 
aspects have been revealed by extensive exploration [T3l HU [15J, [16], [TT], [TBI EE] EUJ 
I2T1 [221 [231 EU [25]. While in two dimensions (2D), a complete understanding of the 
system is still being anticipated. On triangular lattice, many methods have revealed 




(1) 
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that the region 7r/4 < <p < 7r/2 exhibits an AFQ order [2TJ| ¥2T\ |2"8] . On the honeycomb 
lattice, a tensor renormalization group method showed the AFQ order in the region 
7r/4 < (j) < 7r/2 is destroyed by pure quantum fluctuations and there is a transition 
from the plaquette order to the AFM order [29]. A recent work proposed that a three- 
sublattice order exists in the SU(3) point </>o = 7r/4 [30] ■ The quantum Monte Carlo 
simulation [31] found the AFM phase is stable in the region —n/2 < <fi < on a square 
lattice. While in the region of < <fi < 7r/4, there lacks evidence on whether the AFM 
order can survive or not. 

In this paper, we show that a novel type of topological spin liquid might exist on 
a square lattice in the region < < 7r/4, in which a topological edge state circulates 
around the boundary of the system. Time-reversal symmetry is broken spontaneously 
for this non-trivial gapped spin liquid. In a fermion representation we found that the 2D 
spin liquid can be described very well by the (projected) spinless p-wave pairing Bardeen- 
Cooper-Schrieffer (BCS)-type Hamiltonian. We obtained the gapless dispersion of the 
edge modes. For a nonzero longitudinal momentum k x ^ 0, one edge mode splits into 
two half modes that can exist individually, which could be termed £(left) and 7£(right) 
chiral modes, respectively. While for k x = 0, zero edge modes emerge, which turn out 
to be Majorana fermion states [321 [33] . 

The paper is organized as follows. In Section II a brief introduction of the fermion 
representation is presented. In Section III the mean-field theory is introduced and the 
solutions of a gapped chiral spin liquid are obtained. Then the corresponding edge state 
of a ribbon system is explored in Section IV. And a summary is made in Section V. 



2. Fermion representation with hard-core constraint 

Firstly we introduce a fermion representation for quantum spin 1. Each spin has three 
eigenstates \rrij) of with the eigenvalues rrij = — 1, 0, +1. We introduce three fermionic 
operators to generate three independent states (i = a/— T), 

=-L(|m i = -l> + |m i = l», (2) 

fU |0> = Imi = 0) , (3) 

fU\0) =-J=(K = -l)-K = l». (4) 
In terms of / operators, the spin operators can be expressed as 

S? = Kft 2 fi,i-ftifi,2), (5) 

Sf \(hJ).,-.f;J,,). (6) 

s; = Kftih* - ftthi)- (7) 

To restore the Hilbert space of spin 1, the hard-core constraint at each site must be 
imposed, 

3 

E4^ = i- (8) 
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In this way the Hamiltonian Eq. ([I]) is mapped to a frustrated SU(3) fermion model 

m 

H = - J i E : *% F * ■ - J 2 E B li B # + E A * (E f&*# - 1 ) > ( g ) 

(ij) i \ ij. / 

where Ji = cos0 > 0, J2 = cos0 — sin0 > 0, :: denotes normal ordering of operators, 
Aj are the Lagrangian multipliers, and the bond operators are defined as 

Fji = fj t pfi,ni = fj,nfi,i^- (10) 

At=l,2,3 jU=l,2,3 

A similar fermion representation could be found in a recent work |35j . 

3. Gapped spin liquid with time-reversal symmetry breaking 

3.1. Bond-operator mean-field theory 

To find the ground state properties of this spin 1 system, we take the mean field 
approximation by introducing two order parameters for the bond operators Bji and 
F 

(Bi+x,i) — (Fi+y,i) — F, (11) 

-i (B i+x J =B X = Be 1 "*, -i (B i+y>i ) = B y = Be**, (12) 

where F, B, rj x , and rj y are real and to be determined self-consistently. Here we have 
taken a unform phase factor in F, which turns out to be negligible when the mean-field 
equations are established. Two phase factors are kept for B field, and we will see the 
final results only rely on the phase difference At] = % — i] x . Under these prescription, 
the effective Hamiltonian reads 

= E x fU^ + E (toAi. + h-c) + Yl (fl p jifl + h -'-) 

i,ti (ij),n 

- XN A + 2N A J l F 2 + 2N A J 2 B 2 , (13) 

where = — 2J X (Fji), = —2J 2 (Bji), and N A is the total number of lattice sites. 
The hard-core constraint shall be imposed on an average level by minimizing the free 
energy. After performing the Fourier transformation, we arrive at a complex p-wave-like 
pairing of independent flavor of fermions [36J , 

H *ft = \ E ^(k)M(k)^(k) + e , (14) 

k,/i 

so = ^XN A + 2N A J 1 F 2 + 2N A J 2 B 2 . (15) 

where the sum of momentum k is carried out in the first Brillouin zone (1 st BZ), the 
spinor ^(k) = (^fl^, f-k,^, the 2x2 Hermitian matrix 

M(k) = d(k) ■ a (16) 
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with the Pauli matrices a = (a x , a y , a z ) and 

d(k) = (4(k),^(k),4(k)), 
d x (k) = 2J 2 B(cosrj x smk x + cosrj y smk y ), 
d y (k) = 2J2B(smi] x sin k x + sin^ sin k y ), 
d z (k) = A — 1J\F (cos k x + cos k y ) . 

The Lagrangian multiplier Aj is taken to be site- independent, \ 
regarded as a mean field. N\ is the total number of lattice sites, and in fact we have 
defined in this way a mean field Hamiltonian that is similar to the 2D Kiteav model for 
the p + ip superconductors for spinless fermions |33j EZ] • But the hard core constraint 
in Eq. (jBJ) may make it different. By performing the Bogoliubov transformation, one 
can diagonalize the Hamiltonian as 



(17) 
(18) 
(19) 
(20) 

A, which can also be 



H, 



eff 



2 ^ 



*t(k)P(k)^(k)+£ , 



(21) 



where P(k) = u(k)a z , ^(k) = (j^, l-u^j with j Kfl = w k /k,^ - ^k/l kiAt , where the 
coefficients satisfy 



1 + 



2w k fk 



4(k) 
w(k) 
4(k) 
w(k) 
4(k) + id^(k) 



w(k) 

If one chooses a real and even v k (i£ 
The spectrum is 



W (k) = |d(k)| = */^(k)+cE(k) + ^(k). 



(22) 
(23) 
(24) 

f k ), then w k is complex and odd, and vice versa. 

(25) 



The free energy can be worked out as 
3 



P 



E ln d + 



/3w(k)> 



+ #0, 



where /3 



and the ground state energy is 



So = -;^(k) + e, 



(26) 



(27) 



By optimizing the free energy with respect to the mean fields, one obtains the mean-field 
equations as the follows, 
1 
3 



d z (k) Bu(k) 
tanh 



3 1 x ^ 

k 



d z (k) (cos k x + cos k 



(28) 



2w(k) 



"tanh^, (29) 
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B (cosr] x + cos rjy) 



3 1 

2iV! 




d x (k) (sin k x + sin k y ) 



tanh 



(30) 



B (sin^ + sin^) 



3 1 

2 AT 



w(k) 

dj,(k) (sin fc x . + sin fc^) 



tanh 



/Mk) 

2 



2 



(31) 



w(k) 



All the mean fields can be determined by solving the set of mean-field equations self- 
consistently. It is remarkable that all dispersions for quasi-particles are three-fold 
degenerate at the mean field level. Coexistence of non-zero solutions for both mean fields 
F and B affirms the meaningful bond operator decomposition scheme in Eq. (19|)-( TT0|) . 
Notice that the spectrum Eq. (1251) is a function of the phase difference \Ar]\ = \r) y —r) x \, 
we can also take \Arj\ as the optimizing parameter. It is clear to see that the effective 
Hamiltonian Eq. (Tl4l) preserves the time- reversal symmetry when \Arj\ = [34], and 
does not when \Ar)\ ^ 0. The equations are solved at zero temperature to reveal the 
ground state properties. The non-zero solution of mean fields A, F and B for several 
choices of phase difference \Arj\ are illustrated in Figure 1. At zero temperature, one 
can get a simple form for the ground state energy 



where A is cancelled due to the substitution of the mean-field equations Eq. (I28|) -( l3~li) 
in Eq. (1271) . The numerical solution for the ground state energy E and the gap of 
energy spectrum, A gap = min (a; (k)), are illustrated in Figure 2. We found that the 
lowest energy state can be reached by choosing the phase difference |A?y| = tt/2. This 
solution is a p-wave paired gapped spin liquid and breaks the time-reversal symmetry. 
The gapped spin liquid revealed here can be classified by the Z 2 invariant [3E]- And 
by the Z 2 invariant in the topological spin liquid state, we find a special fermion parity 
pattern at high symmetry points in momentum space: even fermion parity at k = (tt, ir), 
k = (0, 7r), and k = (7r,0) and odd fermion parity at k = (0,0) (please see detailed 
calculations in Appendix). The physical regime of the spin liquid should not exceed 
the Heisenberg point 4>h = 0, because the regime with — 7r/2 < < <p H = exhibits 
an antiferromagnetic order [31]. The meaningful numerical solutions with nonzero F 
and B for |A^| = 7r/2 ceases near (p > —0.08 (see Figure 2), a little less than 4>h = 0. 
The discrepancy can be ascribed to the crudeness of the mean-field theory. At the 
SU(3) point 0o = 7r /4 [30], our result shows that the gap closes along the loop line 
k* = (k x , k*) : cos k* + cos k* = 1/A with A = 2.412513447, where the spectrum behaves 
linearly as u (k) ~ c(k*) |k — k*| with anisotropic velocity c(k*). The gaplessness on 
the loop line does not imply any order, so we get a gapless spin liquid that is highly 
degenerate in thermodynamic limit. 

3.2. Ground state and the Chern number 
The p-wave paired ground state of the bulk system reads [36] 




E = —2N\JiF 2 - 2N A J 2 B 2 , 



(32) 





(33) 
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Figure 1. (Color online) Numerical solution for the Lagrangian multiplier A and 
the mean fields F and B. The lines for several selected phase differences are: 
|Af?| = \Wx — Vy \ = 71 black solid line; tt/3, blue dashed line; ir/6, red dotted line; 0, 
green dot-dashed line. Please see more details in the text. 

where the prime on the product indicates that each distinct pair (k, — k) is to be taken 
once. This ground state exhibits a non-trivial topological property that can be signified 
by the Chern number of the spinless SU(3) fermions. For each flavor of spinless fermions, 
the Chern number is defined by [36J EH] 

c = T \ d2k[n (k) ' dk * n (k) x dh * n (k)] (34) 

where D, means the volume of the first Brillouin zone, n (k) is defined as n (k) = nfnxi ■ 
By substituting Eq. (JTZJ) in Eq. (El), we get 

_ sin (At/) / 2J 2 B\ 2 f 2 [^y^ (cos k x + cos k y ) - cos k x cos k y ] 

And by substituting the numerical mean-field solutions at zero temperature in, we obtain 
the simplified result, 

C = ±1, (36) 

for < < 7r/4. Thus the total Chern number of this topological state is C = ±3 due 
to symmetry for the fermions of different flavors. The bulk system's nontrivial ground 
state can be labeled by this Chern number. 
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Figure 2. (Color online) Upper: the ground state energy Eq; lower: the gap of energy 
spectrum, A gap = min(w (k)). The line types are the same as the ones described in 
Figure 1. 



4. Spin edge states 

4-1. Edge modes of the ribbon system 

To demonstrate the spin edge states explicitly, we set up a ribbon (or ladders) system 
with a pair of open edges in y direction and keep periodic boundary condition along 
x axis (Figure 3). Noticing that k x is still a good quantum number, we start with an 
Hamiltonian with L max legs, 

H 'eff = E $ 1 (*») M (*«) ^ (*«) + e , (37) 

k x >0,fi 

®]i(kx) — (/(fc^l),^ —J /(fe,L max ) )A t' /(-fc a ,l),M> f(-k x ,L max ),ti) , (38) 

where we have parsed the zero momentum states in the first term and restricted the sum 
to the positive values of momentum. The matrix M (k x ) are too large to be presented 
here. One can easily solve the Hamiltonian numerically. The resulting diagonalized 
Hamiltonian could be written in the form 

H 'eff = E *m (*«) + (39) 

kx>0,/J, 

P (k x ) = dia [u (kxA) , -,uJ(k x ,L max y, -^(fe^max)] » > W) 
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right 




7 = 1 



Figure 3. (Color online) The ribbon system with open boundaries in y direction. The 
edge modes at the two edges convect in opposite directions. 




Figure 4. (Color online) Numerical solutions of the spectra of the ribbon for the 
selected model parameter </> = 0.435847. The system's size is L max = 50. The lowest 
thick line denotes the edge spectrum Eq. (j46| . Notice the gauge choice A77 = w/2 is 
taken here, so we have right modes for k x > and left modes for k x < 0. The inset 
shows the uppervalue k T for the edge states as a function of model parameter </>. 



^l(kx) — ^7(^,1),^) "•■>llk x ,L mJa ),ftl 7(-fc*,l),^ 7(-fc*,£max),juJ 5 

or a further simplified one 



(41) 



H 'eff= p( t fc Bl i), M 7(fc«,i),A* + T(-fc s ,i), A1 7(-fe,i),/*J- J] w (^)+ £ o-(42) 

k x >0,fi,i k x >0,fj,,i 

We choose the canonical operator 1 - ) (with subscript z = 1) to denote the 
edge excitations above the ground state and ^(^,1) the edge excitation energy. The rest 
modes are bulk modes. The new ground state \Q' Q ) is quite different from that in Eq. 
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(13"3"|) . \Q' Q ) itself contains the chirality of the edges and should satisfy the condition 

7(iW),/u W = °» Vfc ^' A*- ( 43 ) 
We will specify the ground state numerically later when evaluating some quantities. The 
lowest energy mode for each k x could be collected as 

H 'lowest= Y W (fc.,l) 7(^,1)^7(^,1),^ + l\- kxt i)^l{-k x ,l)^ ■ (44) 
k x >0,fj. 

It seems that this effective Hamiltonian denotes the edge states. However, this is not 
necessarily the case. The lowest spectrum is a piecewise function 

/ 2J 2 B\smk x \, (\k x \<k T ); 

U(kx,1) ~ { \J[\ - 2JiF(l + cosk x )] 2 + UlB 2 sin 2 k x , (\k x \ > k T ), ^ 

where k T = axccos(X/2J 1 F — 1) denotes a transition point. At this transition point, 
both co(k x ,i) and its first-order derivative are continuous, but its second-order derivative 
is discontinuous. One would find that only the first piece is the edge mode spectrum, 
i.e. the edge mode spectrum ceases at k T and reads 

w (f ,i) = 2J * B l sin k *\ ' (~ kT ^ k *^ F )' ( 46 ) 
While the second piece is still a bulk mode spectrum. These conclusions are testified by 
our numerical solutions. 

4-2. Numerical solution of the chiral edge states and zero mode Major ana jermion 
states 

Now we discuss some more details about our numerical results. From the numerical 
results, we confirm that the edge modes are well localized at the edges. In practice, the 
solutions with Ar] = tt/2 and Arj = —n/2 are degenerate, but have opposite chirality. 
So we only demonstrate the solution for Arj = it/2 and <p — 0.435847 in Figure 4. For 
< k x < k T , we find the edge mode j}_ kx u ^ localized at the left edge and ^ at 
the right edge. Thus we may call them left (£) and right (1Z) chiral modes respectively, 

C» {~k x ) = 7(-fe,i),n = Y [ U (-k x j)f(k x ,j),v + v (-k x ,j)fl_ kxd) ^ , (47) 

i=i 

^max 

T^n (k x ) = l(k x ,i), M = Y [ U (k x ,j)f(k x ,j),Li + v {k*,j)f\-k x ,j),n ' (48) 

3=1 

where the coefficients U and V are real (as is contrast to the coefficients and for 
the periodic boundary in previous section) and depicted in Figure 5(a) and (b). With 
the increasing size of the system, the zero-momentum edge excitation cof ^ goes to zero 
rapidly. Beyond L max = 50, its energy value is so small that one can hardly decern it 
from the machine precision. We may denote the zero modes as 

^max 

£ M = 7(0,1),/, = Y [ U (oj)f(oj),n + v (o,j)f( ,j),v > (49) 
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where the real coefficients are depicted in Figure 5(c). We see this zero modes manifest 
itself at both edges. It is in fact a Majorana fermion state [32], since the mode contributes 
zero energy to the system and can be rewritten as 

= 1 - 2iMf l M^ (50) 



with two Majorana fermions 



K = E tv,) i Vjm (-'/(o,a, + • < 51 > 

3=1 

imax j r I t r 

= E ^±^> + , (52) 

3=1 

localizing at the two opposite edges. All of the above coefficients satisfy the relations 
numerically for large enough L max (k x > 0) 

El^+^.-J =L ( 53 ) 

3=1 

^(fej) = ^(fc.JJ^C-fc.j) = _V (-W)> ( 54 ) 

U(-k x ,j) = - ^(fc^Lmax-j + l); ^(-fej) = ^(-fe.Lmax-i + l)^ 55 ) 

u (-k x ,j)U(k x ,j) = V { - kx:j) V {kxJ) = 0. (56) 

But notice each mode possess a U(l) symmetry, so that the values may be changed 
according to the symmetry transformation. 

4-3. Edge spin correlation functions 

Although the edge states of the ribbon are clear to see in the fermion representation, 
it is still illusive from the point of view of the spin language. In order to show the 
properties of the spin edge state, we measure the spin correlations and thermodynamic 
quantities, such as the specific heat, contributed by the edge. 

It is well-known a gapped spin system exhibits an exponentially decaying spin 
correlation in the bulk. Of all the spin correlations for the ribbon system, the one at 
the edge is of our great interest. We choose the right edge of the ribbon (Figure 3) to 
measure the spin correlations in the ground state, 

C edge(h % + r ) = C e % e ( r ) = ( 5, fi,L ma x)' S '(i+r,L m a X )) ' ( 57 ) 

Now we need to find out the the ground state \ ft' ). The edge modes in H' edge could be 
singled out and serve as a quasi-lD effective Hamiltonian, 

Heff = Hbulk + Hedges (58) 

Kdge = u ltl) [T(fc X) i) )M 7(fc*,i)^ + l\-k x ,i),un{-k x ,ily\ > ( 59 ) 

0<k x <k T ,ti 

which could be utilized to evaluate quantities along the edges. Near k x ~ 0, the edge 
modes behave linearly wf*^ ~ 2J 2 B \ k x \ [33] (please see the lowest thick line in Figure 
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Figure 5. Visulization of the edge modes in fermion representation on a ribbon with 
width £ max = 50 (1 < j < L max ): (a) the left modes at k x — — n/5, (b) the right 
mode at k x = 7r/5 and (c) the zero momentum mode at k x = 0. The model parameter 
= 0.435847 and gauge choice Ai] — ir/2. 



4). Since the bulk and edge modes are independent, one can write the ground state in 
a separable form 

M)=\tf , bulk )®p' , edge ), (60) 

where \^' 0t buik) an d \^'o,edge) are ^ ne lowest energy states of H' bulk and H' edge respectively. 
We have 

|^0,edge) = Y\. 

0<k x <k T ,j,fi 

where the prime on the product indicates that each distinct pair (k x , —k x ) is to be 
taken once. One can easily verify that 7(^,1),^ \ ^o, e dge) = 0- Since we are concerning 
the quantities along the edge, the edge correlation function Eq. fl57|) can be evaluated 



(kxj) V{k x ,j)J (k x ,j),fiJ (-k x ,j),fi 



J (k-.ft.u.J (- 



10} 



(61) 
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approximately just by \Q' 0edge )- At first, one can work out 



J - 



t 



-U, 



t 



f(k x ,l 



V 2 

{kx max ) 



(62) 



Then the correlation function is deduced as 



- E «■ 



/ /t 



+ 2 



/ 



j-^max) ( )^max)i/^ 



0<k x <k T 



e " (/(l I ,t miI ),^(t,im» I ),/i) 



0<k x <k T 



i-^niax) ^(^a; j-^raax) 



+ 2 



iv A e »w 



0<fco;<fc T 



(64) 



The results show a general power law 

C edge( r ) ~ ^J- 

In practice, we fit the numerical data by the formula 

In C z e l ge (r) rj In a - 5 In r (65) 

instead. At the model parameter = 0.435847, we obtain a ~ 0.0156407 and 
5 ~ 2.174343 (see Figure 6(a)). For other model parameters, the results are not much 
different (Figure 6(b) and (c)). So we see that the edge spin correlations decay like a 
power law along the edge and exclude the possibility of exponential decay behavior in 
the main region of model parameters. 



4-4- Specific heat contributed by the edge modes 

Now we turn to the specific heat. At low temperatures, the bulk states contribute little 
to the specific heat due to the existence of bulk gap, while the gapless edge modes give 
the main contribution. By adopting Eq. (1581) as the effective Hamiltonian, one can 
worked out the contribution of the edge modes to the specific heat that behaves linearly 
in temperature T, 

"Em / m \ 2 



V 



E 



, , , cosh" 2 ( ) p (E) dE « ^-T, (66) 

N A k B J \2k B T ) \2k B TJ HK ! 8J 2 B v ; 

where we have released the upper limit of the integral for simplicity (E m — > oo) and the 
density of states is 



2ttJ 2 B 



(67) 
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0.0 0.1 0.2 0.3 0.4 0.5 0.0 0.1 0.2 0.3 0.4 0.5 



Figure 6. (Color online) (a) The spin correlation functions (in logarithm) along the 
right edge of the ribbon. The dots are the numerical results. Here the width of the 
system is i max = 50, the model parameter is <p = 0.435847 and the gauge choice is 
Ar) = tt/2. The solid blue line is a power-law fit (a « 0.0156407 and S « 2.174343). 
The red dashed straight line is an exponential law fit. (b) The fitted coefficient a as a 
function of cj). (c) The fitted exponent 6 as a function of 5. Please see more details in 
the text. 

5. Summary 

In a brief summary, we showed a possible gapped chiral spin liquid in the 2D square 
bilinear-biquadratic system in the region of < <fi < it /A. As a consequence, the time- 
reversal symmetry breaks spontaneously and an interesting topological ground state is 
revealed. We numerically analysed the resulting spin edge states for a ribbon system 
in detail. This method may be applied to other relating systems to specify a spin 
liquid state. We found £(left) and 7£(right) chiral edge modes for nonzero longitudinal 
momentum k x ^ and a zero edge modes for k x — 0. The power-law decay of the edge 
spin correlation function and the contribution of the nontrivial spin edge state to the 
specific heat at low temperatures are found. In the future work, the properties of low 
energy excitations would be of great interest. 



Topological edge states in spin 1 bilinear-biquadratic model 



15 



Acknowledgement 

The authors thanks Professor Shun-Qing Shen for fruitful discussions. This work was 
supported by SRF for ROCS SEM (20111139-10-2), the Chinese National Natural 
Science Foundation under Grant No.: 11074177, 11174035. This research was supported 
in part by the Project of Knowledge Innovation Program (PKIP) of Chinese Academy 
of Sciences, Grant No. KJCX2.YW.W10. 



Appendix: Classification topological state by Z 2 topological invariants 

For each flavor of fermions (omit the flavor index \i in Eq. (114p ). the effective 
Hamiltonian is 

^// = ^$ t (k)M(k)$(k), (68) 

k 

with M(k) defined in Eq. ({To]) . From the results in Rcf.[38j, the 2x2 Pauli 
matrices can be divided into two groups - even matrix o z and odd matrices, a x and 
o~ y . Because the coefficients of odd matrices are zero at four high symmetry points of 
square lattice in momentum space, we can only focus on the coefficients of even matrix, 
d z (k) = A — 2JiF (cosk x + cosk y ). The four Z 2 topological invariants are defined as 

G = l-0(4(k)), (69) 

where 0(x) = 1 if x > and Q(x) = if x < 0. 

Hence, for points (0,0), (0, 7r), (tt, 0), (tt,tt), the four Z 2 topological invariants are 
explicitly given by 

C M o, ) =e[A-4J x F], 

0c=(0,7T) = 0[A], 

Cfe=(7r,0) = @[A], 

C k=M = e[\ + 2J 1 F). (70) 
For k =(tt, 7r), k = (0, tt), k = (tt, 0), we have a trivial result as 

Ck=(7r,7r) = Cfc=(0,7r) = Cfc=(vr,0) = 0; (71) 

for k =(0, 0), we have 

Ck= ( o,o) = 6 (A - AJ X F) . (72) 

Thus we identify two distinct topological states: the topological state with trivial 
topological invariants 

Ck=(0,0) — Ck=(7r,7r) = Ck=(0,7r) = Ck=(7r,0) = (73) 

for A > AJiF and the topological state with nontrivial topological invariants 

Ck=(7r,7r) = Ck=(0,7r) = Ck=(7r,0) = 0, Ck=(0,0) = 1 (74) 

for A < 4JiF. And in the topological spin liquid state in this paper, we find a special 
fermion parity pattern: even fermion parity at k = (tt,tt), k = (0,tt), and k = (tt,0) 
and odd fermion parity at k = (0, 0). 
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